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We formulate quantum statistical mechanics of particles 
obeying fractional statistics, including mutual statistics, by 
adopting a state-counting definition. For an ideal gas, the 
most probable occupation-number distribution interpolates 
between bosons and fermions, and respects a generalized ex- 
clusion principle except for bosons. Anyons in strong mag- 
netic field at low temperatures constitute such a physical sys- 
tem. Applications to the thermodynamic properties of quasi- 
particle excitations in the Laughlin quantum Hall fluid are 
discussed. 

PACS numbers: 05.30.-d, 05.70.Ce, 71.10+x, 73.40.Hm 



Statistics is the distinctive property of a particle (or 
elementary excitation) that plays a fundamental role in 
determining macroscopic or thermodynamic properties 
of a quantum many-body system. For example, super- 
fluidity or superconductivity is essentially due to Bose- 
Einstein condensation; and the stability of macroscopic 
matter PJ depends crucially on the Fermi-Dirac statis- 
tics of electrons and protons. In recent years, it has been 
recognized that particles with "fractional statistics" in- 
termediate between bosons and fermions can exist in two- 
dimensional H or in one-dimensional HQ] systems. Most 
of the study has been done in the context of many-body 
quantum mechanics. Despite calculations of certain ther- 
modynamic properties in a few examples |3|,^p0 with the 
help of exact solutions, general formulation of quantum 
statistical mechanics (QSM) with a distribution that in- 
terpolates between bosons and fermions is still lacking. 

A single-particle quantum state can accommadate an 
arbitrary number of identical bosons, while no two identi- 
cal fermions can occupy one and the same quantum state 
(Pauli's exclusion principle). In QSM 0, this difference 
gives rise to different counting of many-body states, or 
different statistical weight W . For bosons or fermions, 
the number of quantum states of N identical particles 
occupying a group of G states is, respectively, given by 



W b = 



(g + n-i)\ 



or Wf = 



G\ 



Nl (G - 1)! ' 1 N\ (G- 

A simple generalization and interpolation is 
[G+(iV-l)(l-a)]! 



N)\ 



W = 



N\ [G - aN - (1 - a)]\ 



(1) 



(2) 



with a = corresponding to bosons and a = 1 fermions. 
Such an expression can be the starting point of QSM for 
intermediate statistics with < a < 1. Let us first clarify 



its precise meaning in connection with "occupation of 
single-particle states" , and justify it as a new definition 
of quantum statistics, a la Haldane 0]. 

Following ref. Q, we consider the situations in which 
the number Gi of linearly independent single particle (or 
elementary excitation) states of species i, confined to a 
finite region of matter, is finite and extensive, i. e. pro- 
portional to the size of the matter region in which the 
particle exists. Now let us add more particles with the 
boundary conditions and size of the condensed-matter re- 
gion fixed. The A^-particle wave function, when the coor- 
dinates of the N — 1 particles and their species are held 
fixed, can be expanded in a basis of wave functions of the 
rt-th particle. The crucial point is that in the presence of 
other particles, the number di of available single-particle 
states in this basis for the n-th particle of species i gen- 
erally is no longer a constant, as given by Gf, rather it 
may depend on the particle numbers {-/V,} of all species. 
This happens, for example, when localized particle states 
are non- orthogonal; as a result, the number of available 
single-particle states changes as particles are added at 
fixed size and boundary conditions. Haldane |jj defined 
the statistical interactions ctij through the relation 



Ad % = -J2< 



'■ij AiV,- 



(3) 



where {ANj} is a set of allowed changes of the parti- 
cle numbers. In the same spirit, but more directly for 
the purposes of QSM, we prefer to define the statistics, 
including possible mutual statistics, by counting the 
number of many-body states at fixed {A^}: 



W = 



n 



G. ■ A, 1 )J f n,,,':.V„ <*,-,)]! 
i (Ni)\ [G.-l-V.a^-^-)]! 



(4) 



The statistics parameters onj must be rational, in order 
that a thermodynamic limit can be achieved through a 
sequence of systems with different sizes and particle num- 
bers. If we restrict to only one species, ([I]) reduces to (Q). 
Note there is no periodicity in so-defined statistics a. 

Now let us consider the simple case that every single- 
particle state of species i has the same energy e^. An 
ideal gas, by definition, is a system whose total energy is 
a simple sum 



E = J2N iSi 



(5) 



where each term is linear in the particle number A^. De- 
fering the discussions about when eqs. (§ and (§ are 



1 



obeyed, let us first apply them to study QSM of such 
ideal gas. Following the standard procedure ||], one may 
consider a grand canonical ensemble at temperature T 
and with chemical potential pn for species i. According 
to the fundamental principles of QSM, the grand parti- 
tion function is given by (with k the Boltzmann constant) 

Z = E W ({ N *}) exp{^A^ ~ e l )/kT} . (6) 

{#«} * 

As usual, we expect that for very large Gi and N%, the 
summand has a very sharp peak around the set of most- 
probable (or mean ) particle numbers {ATj}. Using the 
Stirling formula log AT! = Alog(A/e), and introducing 
the average "occupation number" defined by rii = Ni/Gi, 
we express logW as (with /3y = ctijGj/Gi) 

^2G l {-n l \ogn t - (1 - ^2 fan j)log(l -^PijTij) 

i 3 3 

+ [1 + ]T(% - 0ij)nj] log[l + - PM)- (?) 

3 3 

The most-probable distribution of n,; is determined by 



^- [log W + 2 G ini {m - s l )/kT] = , 

i 

It follows that 
n . e ( £i -w)/feT = (1 + ^ (5 . fc _ p ik)nk] 



(8) 



n[ r 



1 - Efe ftfc^ 



Setting Wi = n l 1 - (3 ik n k /ni, we have 



(9) 



(10) 



Therefore the most-probable average occupation num- 
bers rii(i = 1, 2, • • •) can be obtained by solving 



(11) 



with uii determined by the functional equations (|10|). The 
thermodynamic potential f2 = — kT log Z is given by 

. „ 1 + m — Y"* ■ /3i»rii 

n^-PV = -kT^G t \og V* 3 ; (12) 



and the entropy, 5 = (J5 - Vi N i ~ ty/ T , ^ 



! = !>{• 



log 



1 + ni - V ( Aj^i 



j A^U 71 3 



}. (13) 



Other thermodynamic functions follow straightforward- 
ly. As usual, one can easily verify that the fluctuations, 



(N^ — Aj 2 )/Aj 2 , of the occupation numbers are negligi- 
ble, which justifies the validity of above approach. 

To see more explicitly the consequences of the above 
results, consider the simplest case of only one species with 
different energy levels, for which we set aij — aSij and 
/ij = /.t. Then the average occupation number rii satisfies 



(1 - am) a [l + (1 - a)rn\ 



me 



and we have the statistical distribution 

1 



(14) 



(15) 



1 w(e(^-^/ kT )+a ' 

where the function w(£) satisfies the functional equation 

w(O a [l + w(0] l - a = ( = e^V kT . (16) 

Note that w(() = C — 1 for a = 0; and w(() — £ for 
a = 1. Thus, eq. (jig ) recovers the familiar Bose and 
Fermi distributions, respectively, with a = and a = 1. 
For semions with a = 1/2, eq. ( |l4]) becomes a quadratic 
equation, which can be easily solved to give 

(17) 



v/1/4 + exp[2( £l - ,i)/kT\ 

For intermediate statistics < a < 1, it is not hard to 
select the solution iu(C) of eq. (0) that interpolates be- 
tween bosonic and fermionic distributions. In particular, 
when £ is very large, we have u>(C) ~ C and, neglecting a 
compared to w(Q, we recover the Boltzmann distribution 



(18) 



at sufficiently low densities for any statistics. 

Furthermore, we note that £ is always non-negative, so 
is w; it follows from eq. (15) that 



rii < l/a 



(19) 



This expresses the generalized exclusion principle for 
fractional statistics. In particular, at absolute zero, £ = 
if Si < (i, and £ = +oo if Si > fi. From eq. (|l6|), we have 
w = and oo respectively. Thus, we see that at T = 0, 
for statistics the average occupation numbers for 

single-particle states with continuous energy spectrum 
obey a step distribution like fcrmions: 



rii = 0, if £; > Ef; 
l/a, if Si < E F . 



(20) 



The Fermi surface e = Ep is determined by the require- 
ment J2 e <e f ~ a -N ■ Below the Fermi surface, the 
average occupation number is l/a for each single-particle 
state, a consequence of generalized exclusion principle. 

One may be tempted to consider, in parallel to usual 
Bose and Fermi ideal gas, the case Q with 



2 



Si 



2m ' 



G* 



2^ 



(21) 



say in two dimensions with V the area. Then treating 
momentum as continuous, one has the Fermi momentum 



hp 



(4ira)(N/V). 



(22) 



Moreover, at finite temperatures, by using (U) and (fTeJ) 
the sum J2i Gifii = N can be done to give 



JL 

kT 



2irh 2 N 

i 

mkT V 



log 



1 — exp 



Using the identity derived by integration by parts, 
1 — an; 



dsi log - 
,, 1 + (1 - a)rii 



dEiE; 



(24) 



we have the statistics-independent relation PV = E. In 
the Boltzmann limit (exp(p/kT) << 1), w(() = £+a — 1, 



PV = NkT[l + {2a - 1)NX 2 /W] , 



(25) 



where A 



2nh 2 /mkT. So the "statistical interac- 



tions" are attractive or repulsive depending on whether 
a < 1/2 or a > 1/2. 

Whether a given system satisfies the seemingly harm- 
less conditions (0) or (^lj) together with (0) is a nontrivial 
question. The state-counting definition (0) for fractional 
statistics does not apply to free anyons, i.e. Newtonian 
particles carrying flux-tubes fll]Jl^] • On the other hand, 
anyons in magnetic field satisfy (0T, but not the condition 
(|J) unless (at very low temperatures when) all anyons are 
in the lowest Landau level (LLL) O] . The Jastrow-type 
prefactor U a< b(z a — zi,) ^ , with < 6 < 2tt, in the anyon 
wave function has the effect of increasing the flux through 
the system by (8/tt)(N — 1). Thus, with fixed size and 
number of flux, the dimension of the effective boson Fock 
space [|l4|l5 is given by d = — (9/n)(N — 1), where 
= qBV/hc = V/Vo, with q the anyon charge. Eq. (0) 
applies, with the single-anyon degeneracy G = and 
the statistics a = 9/tt. Applying eqs. ( |l5|) and ( ^6|) with 
only one energy e = hu> c /2, we have 



N 
G 



P_ 

Pa 



; ( e (e-M)/fcT) + , 



(26) 



where u>(C) is the positive solution of eq. (|l6[). Here 
p = N/V is the areal density, and po = 1/Vb- Eq. (26), 
together with (|l6|), determines the chemical potential p 
in terms of p/po and T. Thermodynamic quantities can 
also be expressed as functions of the ratio p/po- In par- 
ticular, the thermodynamic potential is 



V a w 



-*r£lo 6 I±<I^. ,27) 
Vq 1 — an 



PV 
NkT 



^Po J 



log 



l + (l-q)(p/p ) 
1 - a(p/p ) 



(28) 



The pressure P is linear in T for fixed p. It diverges at 
the critical density p c = (l/a)po, which corresponds to 
the complete filling of the LLL. The emergence of an in- 
compressible state at filling fraction 1/a is a consequence 
of the generalized exclusion principle ( pp| ) . Jl6| The mag- 
netization per unit area is 

2u , 1 + (1 — a)(p/p Q ) . N 

M = -p p+^io g — ± — , ( 29 ) 

A 1 - a(p/po) 

where /xo = qh/2mc is the Bohr magneton. Note the first 
(de Haas- van Alphen) term is statistics- independent. At 
low temperatures, kT << ?iw c , the second term can be 
neglected except for p very close to (l/a)po, where it 
gives rise to a non-vanishing, a-dependent susceptibility 



X 



kT- 



U 
B' 



p/Po 



2irhc y B' (l-ap/po)[l + (l-a)p/p ] 
The entropy per particle is also a-dependent: 



(30) 



S_ 

N 



= fc(l — a 

-/clog — 
Po 



Po, 

p' 



log[l + (l-a)-^] 
Po 



k(^ 
P 



a) log(l — a — ) 
Po 



(31) 



Eqs. @, (||) and @ have been derived in Ref. § 
from the known exact many-anyon solutions in the LLL. 

Vortexlike quasiparticle excitations in Laughlin's in- 
compressible 1 /m-fluid (m being odd) JlT[ | are known to 
be fractionally charged anyons, and their wave functions 
are such as if they are in the LLL (with electrons acting 
as quantized sources of "flux") Q- 0. The existence of 
two species of excitations, quasi-holes (labeled by — ) and 
quasi-electrons (labeled by +), dictates nontrivial mutual 
statistics. From the "total-flux" constraint 

ee eBV/hc = mN e + iV_ - N + , (32) 

where N e the number of electrons, it follows Q that 

a = — = a |_ = 1/to, a++ = 2 — I/to , (33) 

(appropriate for hard-core quasielectrons ]l5| , pl"| | ) . The 
single excitation degeneracy in the thermodynamic limit 
is G+ = G_ = (l/m)N c f ) . Ignoring the interaction en- 
ergies and assuming the system is pure, we apply the 
formulas ( |l0(|l5| ) to this case with the values (|3^) for a y - . 
The densities p± of the excitations are given by 



P± 



w T + a TT — a± T 



(34) 



po (w+ + a ++ )(w- + a — ) - a + _a_ + 
where po = G±/V; w± satisfy the functional equations 

W± "±± (1 + W± )!-«±± f ° T± = e (=±-P±)/fcT . 

V 1 + J 



The equation of state is 



(35) 



3 



Here e± is the creation energy of a single excitation. At 
T = or very close to it, thermal activation is negligi- 
ble; there is only one species of excitations, behaving like 
anyons in the LLL. Say, if < mN e , then there are 
only quasielectrons: n_ = 0, n + = (mN e — NA/G+ = 
l/(u)+ + One can apply the above eqs 



At higher temperatures, thermal activation of quasipar- 
ticle pairs, satisfying + /i_ = 0, becomes impor- 
tant and the effects of mutual statistics become manifest 
with increasing density of activated pairs. The thermo- 
dynamic properties at different sides of electron filling 
v = NpjN,^ = 1/m are not symmetric due to asymme- 
try in quasielectrons (a++ = 2 — 1/m) and quasiholes 
(a = 1/m). The general equation of state is 



— = Po log 



kT 



1 + {.Pi/ Pa) - Ej a i](Pj/Po) 
1 - E, a ij(pj/po) 



When the excitation densities satisfy 



E 

]=+,- 



iPj = Po 



(36) 



(37) 



the pressure diverges and a new incompressible state is 
formed, as a result of the generalized exclusion princi- 
ple obeyed by the excitations upon completely filling the 
LLL. At T = 0, e.g., for i = + (quasiparticles) the limit 
( |37| ) is reached when electron filling is v = 2 /(2m — 1), 
giving rise to the well-known hierarchial state |l4|,|22|,|l5| . 
At finite T, it may happen at somewhat different filling, 
because of the additional quasihole contribution in eq. 
( |37| ) . Moreover, the magnetization per unit area is 



-PiPi 



ekT , po 
— — log — 

mnc 



+ Pi - Ej Q-jjPi 

PO - E 7 OijPj 



(38) 



Here //± = de±/dB; the linearity of e± in B is observed 
in recent experiments |23| . Hopefully, when kT is of or- 
der of e± or higher, the ay-dependent second term may 
give an appreciable contribution. The mutual-statistics- 
dependent total entropy, S = Ej is 



Si 
k 



log- 



Po 



Po 



E 



am — )log(l - y^Oj 



ft 



ft- 



+ [i + - E f, '<- ^st 1 + E(*« - ( 39 ) 

Pi p. Po 

To conclude, we have formulated the QSM of particles 
of fractional statistics (and mutual statistics), adopting 
a state-counting definition. When their interaction ener- 
gies can be neglected, they obey a statistical distribution 
that interpolates between bosons and fermions, and re- 
spect a generalized exclusion principle which makes them 
behave more or less like fermions. Anyonic quasiparticle 
excitations in the Laughlin quantum Hall liquids are such 
particles with nontrivia! mutual statistics. It would be 



interesting to experimentally test their thermodynamic 
properties. Also it would be interesting to see if and how 
statistics with a > 2, which makes perfect sense with the 
state-counting defintion (^), realizes in Nature. 
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